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A framework is developed to describe the two-point statistics of potential vorticity in ro- 
tating and stratified turbulence as described by the Boussinesq equations. The Karman- 
Howarth equation for the dynamics of the two-point correlation function of potential vor- 
ticity reveals the possibility of inertial-range dynamics in certain regimes in the Rossby, 
Froude, Prandtl and Reynolds number parameters. For the case of large Rossby and 
Froude numbers, and for the case of quasi-geostrophic dynamics, a linear scaling law with 
2 /3 prefactor is derived for the third-order mixed correlation between potential vorticity 
and velocity, a result that is analogous to the Kolmogorov 4/5-law for the third-order 
velocity structure function in turbulence theory. 



1. Introduction 

Potential vorticity q = u> a ■ Vp, where u: a is the total v orticity and p is the density, 
is a lagrangian invariant in rotating and stratified flows fertell |l942)). In the quasi- 
geostrophic (QG) limit, in which the rotation and stratification effects are strong, the 
potential vorticity evolution entirely determines the dyna mics of quantiti e s such as wind 
speed , pressure and te mperature fields (se e for e xample. iHoskins et oil (^985); R hinesl 
l|l98r1) : lMulierl (|199!t|L - lHavnes fc Mclntvrel l(l99(tl ). Instead of studying q as a local in- 
variant in the traditional way, we here focus on the global statistical description of the 
potential vorticity field in the s tyle of the Kolmog orov 1941 theory for the statistical 
description of the velocity field ijKolmogorovl l)l94l|) ). We treat potential vorticity as a 
statistical variable and consider the dynamics of its second-order moment in the triply 
periodic or infinite domain. Thus the potential enstrophy Q = (q 2 )/2 is the conserved 
quantity of interest, where (•) denotes a suitable average (over ensembles or over the flow 
d omain) . 

ICharnevI l|l97lj) showed analytically that in the QG limit, conservation of energy and 
potential enstrophy implies an inverse cascade of energy with energy spectrum scaling 
as E(k) ~ fc -5 / 3 for the low wavenumbers, and a forward cascade of potential enstrophy 
with E(k) ~ fc -3 at high wavenumbers. The strong constraint on the energy by Q is 
analogous to the constraint on energy by enstrophy in 2D turbulence, resulting in the well- 
known inverse cascade of energy, the forward cascade of enstrophy, and the corresponding 
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scalings of the energy spectrum l|Kraichnanl 1119711^ . iHerring e.t. nil ill 99^ performed 
moderate-resolution, spectral simulations of non-rotating, non-stratified turbulent flow 
with a passive scalar 6, in which the Ertel potential vorticity q = u> ■ V6> with uj = V x u. 
They studied the evolution of the potential enstrophy spectrum and demonstrated that 
the transfer terms and the dissipation terms are not separated in wavenumber space. 
They argued therefore, that an 'inertial' range of scales dominated by the flux of potential 
enstrophy does not seem possible. The concept of an inertial range is a cornerstone of 
the K41 theory of turbulence. Therefore it seemed futile, in a certain sense, to further 
explore statistical approaches for potential vorticity dynamics. 

To revisit the possibility of potential enstrophy intertial-range dynamics, we depart 
from previous approaches in two ways. First, we confine ourselves to statistics in physical 
space, and undertake a novel study of the evolution equations for the two-point spatial 
correlation function of potential vorticity. Second, we examine, in addition to the QG 
and non-rotating, non-stratified limits, four other regimes of interest in the rotation 
and stratification parameters. Our approach has revealed that some of these regimes do 
allow for the existence of an inertial range of potential enstrophy. Furthermore, when 
we take a prescribed sequence of limits in the non-dimensional parameters and assume 
local isotropy, then exact scaling laws for small-scale potential vorticity statistics emerge. 
These new laws are analogous to existing laws that form the backbone of statistical 
hydrodynamics in 3D non-rotating, no n-stratified flows. 

Following K41 turbulence theory, the lvon Karman fe Howartbl l|l938j) equations for the 
two-point correlation functions of velocity lead to the concept of an 'inertial range' of 
scales dominated by net downscale flux of kinetic energy by nonlinear transfer with little 
dissipative loss. The inert ial-range concept was extended to th e flux of pa s sive s calar 
vari ance bvlYagloml lll949[). and more recentl y to helici t y flux bv lChkhetianil 1 199(t) (see 
also iL'vov e.t all l|l997h : IComez et all (j200r^ : iKurienl 1)20031) V Under the assumptions 
of statistical homogeneity and statistical isotropy of the small scales (local isotropy) , all 
these results take the form of exact scaling laws for the appropriate two-point third-order 
correlations as follows: 

((u L (x + r) -u L (x)) 3 ) = —er; e = 2v(SJu ■ Vu) (1.1a) 

5 

({u L (x + r) - u L (x)){8{x + r) - 6(x)) 2 ) = -^e e r; eg = 2a(V6 ■ V8) (1.16) 

2 

{(ul(x + r) — ul(x))(ut(x + r) X Ut(x))) = — hr 2 ; ft = 2v(Vm • Vw) (1-lc) 

15 

where r is the separation vector between the two measurement points; the velocity com- 
ponent along r is called the longitudinal velocity ul = u r (we will also refer to its vector 
representation ul — u^r); the component orthogonal to r is called the transverse veloc- 
ity ut = u—ul t; 9 is the passive scalar. The quantities e, eg and h are, respectively, the 
mean dissipation rates of kinetic energy, passive scalar variance and helicity, and are de- 
fined in terms of the viscosity v and the thermal diffusion coefficient a. The kinetic energy 
and passive scalar variance have positive-definite dissipation rates, while helici ty does not. 
The la t ter fact, however, do es not in itself preclude a helicity inertial-range l|Chen et all 
( 2003 UKurien et all l|2004j) 'l. Equations (JTTJ are three of the few exact, nontrivial re- 
sults known in the theory of statistical hydrodynamics. They are valuable b enchmarks for 
the st u dy of high Reynolds number ( Re) turbulence in both expe riments l|Dhruva et all 
lll997h : lMvdlarski fc Warhaftl Jl998tk IChambers fc Antonial lll984h\ and numerical simu- 
lations l(Sreenivasan et all lJl99fiMCTOtoh et all (|2002h : iTavlor et all l|2003l) : lKurien et all 
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In section|21we derive the Karman-Howarth equation for the two-point correlation func- 
tio n of q starting from t he Boussinesq equations for rotating and stratified flows. As noted 



bv lHerring et all ( 199^) for zero rotation and stratification, the flux and viscous-diffusion 
terms are not necessarily separated in scale (or wavenumber). However, in statistically 
homogeneous flow with Re — > oo and fixed Prandtl number Pr, there exists the pos- 
sibility of a range of scales wherein the viscous-diffusion rate is at least sub-dominant, 
if not negligible, compared to the transfer rate. Assuming such a range of scales exists 
at sufficiently high Re, we recover a balance between the divergence of the third-order 
correlation between velocity and potential vorticity, and the production-dissipation rate 
of potential cnstrophy. 

In section|3we consider different limits of the Rossby (Ro) and Froude {Fr) numbers. 
When rotation and stratification are small [Ro and Fr are large), we can reasonably 
hypothesize local isotropy at small-enough scales. Then, an exact scaling law analogous 
to Eqs. <|1.I[) arises for the third-order mixed correlation of u and q, with a prefactor 
of 2/3. In the QG limit {Ro and Fr are small), the dependence on Ro and Fr may be 
scaled out in a coordinate system with vertical stretching. For QG scales in the stretched 
coordinates, we again assume local isotropy and deduce the same linear scaling law for 
the mixed third-order correlation. In section^ we tabulate several intermediate limits in 
Ro and Fr; in some cases the correlation dynamics allow for a scale separation between 
flux terms and viscous-diffusion terms, and hence the clear possibility for an inertial 
range of scales dominated by potential enstrophy flux. 



2. The Karman-Howarth equation for potential vorticity q 

The Boussinesq equations for rotating, stably stratified and incompressible flow are: 

—u + fzxu+—X7p+-?-gz = v\7 2 u, V-ii = 0, —p-bw = nV 2 p, (2.1) 
Dt po p Dt 

where D/Dt = dt + {u- V), u is the velocity, w is its vertical component, p is an effective 
pressure, / = 2f2 is the Coriolis parameter, Q is the constant background rotation rate, 
v = p/po is the kinematic viscosity and n is the mass diffusivity coefficient. The total 
density is p{x) = po — bz + p(x), where po is the constant background, b is also constant 
and larger than zero for stable stratification, and p is the density fluctuation such that 
|p| < \bz\ <p - 

We begin with the equation for q in rotating, stably stratified flow l|Embid fc Maidal 
l)l998|) 1 with periodic boundary conditions: 



where u> = us + fz is the absolute vorticity and u> = V x u is the relative vorticity. The 
equation for the potential enstrophy Q = q 2 /2 is obtained by multiplying Eq. 12.21) by q. 
The mean production-dissipation rate of Q over the domain is then given by: 




D 1 D / VT \ 

— — — = (0J a ■ VP) 

Dt DV HJ 



'a 



(2.2) 



— (Q) = v{q V 2 w ■ Vp) + K (q V 2 S7p ■ u a ) = -eg 



(2.3) 



where (•} denotes volume integration over the periodic domain. The quantity sq is not 
sign -definite, allowing fo r both production and dissipation, of potential enstrophy (see 
also lHerring et all (jl994lU . 
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2.1. Equation for the two-point correlation function of q in homogeneous flow 
Following lEmbid fc Maidal l)l998j) . one may define non-dimensional parameters: 



Lf LN 



r = 



BgL 



Pr=~, A = TFr 2 , Re = — (2.4) 



where U, L, L/U are the characteristic velocity, length scale and time scale, and Bpo is 
the characteristic scale for the density fluctuations (B is a dimensionless constant). The 
buoyancy frequency N is given by N — (gb/ po) 1 ^ 2 ■ For the limiting cases Fr — > and 
Fr — > oo, conservation of total energy (u ■ u + p 2 )/2 requires that T = 1/Fr, and thus 
A = Fr. The latter equality is assumed throughout this work. Then the characteristic 
velocity is given by U = Bg/N = B(p g /b) 1 / 2 . 

With the above definitions, the non-dimensional form of the equation for q as measured 
at point x becomes 

^-[u> ■ V/5 + Ro- 1 ^ - Fr^uj 3 } = Re~ l (Fr^V 2 u>z - Vp ■ V 2 uj) 
Dt dz 

-(itePr)- 1 (iZo- 1 V a ^+«-V a V^) (2.5) 
Oz 

where all variables are non-dimensional and uj 3 is the vertical component of u>. From now 
on we include in q only the contributions including the fluctuations oj and p 



q = u-Vp + Ro^ 1 ^ - Fr- 1 ^, 
oz 



(2.6) 



and we work entirely in the non-dimensional units. 

One may write down an equation identical to (|2.5(1 for q', the potential vorticity at 
point x' . Cross- multiplication and summing of the two resulting equations yields the 
equation for the two-point quantity qq' — q(x)q(x'): 



di^^dl 



Re- 



(widip)(wjdj>p') + Ro l {^-w' i d v p' + ^jujidip} - Fr l {u)zu' i d v p' + w' 3 Uidi 



^9~P9~P' R -^Fr-^%+^} + Fr- 2 ^ 



Ro 

Oz Oz' 
(u; k d k p)(d 2 ,ujl)(d t ,p 



Ro 



dz' J <9z 

Ufa 

dp' 



Re~ l Pr 



Fr-^co^Lu'^p') + w'^UiWp) + {Widiflifywti + (^a^')(9> 3 )} 

Ro- l Fr- x &}^ + %d 2 u,} + Fr- 2 {^d 2 ^ + u 3 d]w 3 } 
Oz J Oz J j j 

(uApHided 2 ?) + {^dyp'^iididlp) 



-U d P, Ji 



dp 



Ro-^u'M^p') + ^(did'p) + (uidipM^) + (cjid.p'Kd 2 ^)} 
Oz J oz oz' J oz 

Fr- 1 {cj 3 L>' i (a i ,8§,p') + ^i{did 2 p)} 



Ro 



dz 3 dz' dz' 3 dz 



^di } 



where d{qq')/dt is defined by d(qq')/dt = d t {qq') + q'uidiq + qu'^q 1 . Here all primed 
variables denote their values at x' and cV denotes differentiation with respect to x' v 
Next one may express the equation in terms of the separation vector r between the 
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points x and x', and rewrite the equation in te rms of the new independent variables 
r = x' — x and X = (x + x')/2 (see for example iHil J (Hqq2)), where 

di = -dr i + \dx i ; d v =d ri + ^d Xi ; df = d? = d% + \d Xi . (2.8) 

The procedure to derive t he equation for th e two-point correlation function is familar 
from, for example iFrischl l|l995f) : iHilll l|2002|) . One can perform a change of variables 
and then ensemble average the equation assuming statistical homogeneity. Ensemble 
averaging commutes with the derivative operations d Ti and d Xi ■ Statistical homogeneity, 
as in the case of periodic boundary conditions with constant N and /, implies that the 
derivative operation d Xi acting on any statistic yields zero. The result is 

^f 1 - ~ d rt (q q'( Ui - u{)) - Re- l d ri (qp'd^i - q'pdfa) 

+ {RePry^rM^dlP' ~ faattfp), (2.9) 

where all factors of Ro and Fr are 'hidden' in the non-dimensional expressions for q and 
q' defined by Eq. (|2.6fl . p = —Fr~ 1 z + p and u a = Ro~ x z + u>. Relation 1)2. 9[) is the 
non-dimensional Karman-Howarth equation for the two-point, second-order correlation 
function of q. It describes the dynamics of the two-point correlation of q in a statistically 
homogeneous Boussinesq flow, and is the starting point to consider various limiting cases 
as described below. 



2.1.1. Non-diffusive small-scale limit in homogeneous flow 

For finite Re and Pr, the viscous-diffusion terms on the right-hand side of Eq. 1)2.9)1 
may not, in general, act only at scales much smaller than the non-linear (transfer) terms 
on the left-hand side. Furthermore, each individual term on the right-hand side is not 
sign-definite, and thus allows for both production and dissipation of (q q') at a given 
scale. The latter is reminiscent of the sign-indefiniteness of the dissipation rate of helicity 
and helical velocity statistics; despite this, it has recently been shown that an inertial 
range of helicity can exist with an e x act scaling law recove r ed theoreticall y bv lChkhetianil 
l|l996ft (see als olL'vov et all dl997h : iGomez et all l)2000|) : iKurienl l)2003^ . as well as in 
simulations bv lKurien et all (|2004f) . 

We proceed with a particular order of limits, first Re — » oo with Pr fixed, followed 
by r — > 0. The expectation is that this procedure will access a range of (sufficiently 
small) scales such that the effect of the diffusion is sub-leading compared to the transfer 
for a given scale. One may then derive the leading-order balance between the flux and 
production-dissipation terms. In the first limit, the right-hand side of Eq. I|2.9|l becomes 
negligible for given r. Then, in the limit as r — > 0, d t (qq') — » 2d t {Q) = ~2sq and hence 
Eq. (|2.9jl reduces to 

V r ■ (q q'(u - u')> = -2e Q . (2.10) 

Implicit in this step is the assumption that the production-dissipation rat e ep remains 
finite in the non-diffusive limit Re — * oo, Pr fixed. This is analogous to the iKolmogorovl 
dl94lh hypothesis that kinetic energy dissipation rate remains finite in the inviscid limit. 
Similar assumptions were used for the rates of dissipation of passive scalar variance 
and helicity to derive Eqs. 1)1. ip . Equation (|2.10() is analogous to the result for third- 
order velocity structure functions in stati s tically homogeneous, high Reynolds number 
turbulence in the small scale limit IjFrischl 1)19951^ . 
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3. Limiting cases in the rotation and stratification parameters 

Our starting point is the Karman-Howarth equation for homogeneous flows Eq. (|2.9|1 , 
upon which we impose limits in Ro and Fr using Eq. (j2.6|) (see also Eq. I|2.7|l ). In the 
large Ro, large Fr limit, we derive an exact balance for the isotropic small scales between 
sq and the mixed third-order correlation of u and q. The scaling law thus derived, which 
we will call the '2/3-law', is analogous to the established scaling laws for kinetic energy, 
passive scalar variance and helicity in Eqs. Ql.l[l . The same 2/3-law is shown to hold for 
the QG limit in a stretched coordinate system, presumably in a range of larger scales for 
which vertical stretching leads to local isotropy. 



3.1. Rotation and stratification are small (large Ro, large Fr) 

The equations for large Ro and large Fr should tend toward the equations for incom- 
pressible, variable density 3D flow. Consider the following scalings: 

IN 1 

Ro= — -, Fr=-, Re^O(l), Pr fixed as e -» 0. (3.1) 
e / e 

After the change of variables according to (|2.8|l . ensemble averaging and assuming ho- 
mogeneity, the leading-order 0(1) balance of H2.9J) is 

- d Ti {q q'( Ul - u[)) = Re-^drAq /5'V' 2 ^) - d n {q' fN 2 ^)] 

+ Re- 1 Pr- 1 [d ri (qu' i V 2 p') -d^iq'uj^p)], where q ~ w • Vp. (3.2) 

Note that, unlike in the equation for velocity correlations, the order of the derivative 
with respect to r is the same for the flux and viscous-diffusion terms. Thus we cannot, 
without further assumptions, expect an 'inertial-range' transfer dominated by the flux 
of potential enstrophy in some ra nge of scales. I n this limit, the equations of motion are 
identical to those investigated bv lHerring et all (^994), that is, the density is a passive 
scalar and the momentum obeys the Navier-Stokes equations. 



3.1.1. Local isotropy 

In the large Ro and large Fr limit, it is reasonable to further assume local isotropy 
for sufficiently small scales, that is, invariance of the correlation tensors under arbitrary 
rigid rotations. Invariance with rotation by ir radians about each of the coordinate axes 
yields the constraint that only the longitudinal components of the tensor correlations are 
non-zero. Then Eq. I|3.2|l reduces to 

-V r -(q q'(u L - u' L )) = i?e- 1 V r • (qp'd} u' L - q'p8^ L ) 

+ (Re Pry^r ■ {qu' aL d2p" - q'u; aL d?p), (3.3) 

where subscript L denotes the longitudinal component. Equation (|3.3|l is a Karman- 
Howarth equation for two-point potential vorticity statistics in 3D incompressible, vari- 
able density, statistically homogeneous and locally isotropic flow. We may use the fol- 
lowing isotropic forms for the scalar and tensor correlations in Eq. (|3.3(l : 

r ■ 

(a q') = (q(x)q{x + r)) = C(r) ; (q(x)q(x + r)(ui(x) - m{x + r))) = F{r)^-; 
(qp'd 2 ^ - q'pdfa) = Grirf-i ; (qu/^p? - q'u CH d)p) = G 2 (r)^ (3.4) 
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where C(r), F(r), G\(r) and G2(r) are scalar functions of r. Substituting these isotropic 
forms into Eq. (|3.3(1 and using the identities 

dr Ti _ d dr d Ti d 
dri r ' dri dri dr r dr 

we find 

^ ~ hU^) - «*-^{r>G l[r j) + (a iVl-^pW)). (3.6) 



Relation 1)3.6)1 is the the scalar form of the Karman-Howarth equation for the second- 
order moment {qq'} in locally isotropic flow. 

3.1.2. Statistically steady state, non-diffusive and small-scale limit 

The analogy to the Kolmogorov 4/5-law for longitudinal structure functions may be 
recovered in the non-diffusive limit using the following steps. First, assume a statistically 
steady state in time. Then, in order to observe the limit where the viscosity and mass 
diffusion contributions are small, take the limit Re — > oo with Pr constant, eliminating 
the right-hand side of (|3.6(l . Follow this with r —> to obtain: 



/ dC{r) 



dt K nr ~ v dt \ d 



r 



1 d 



CWIr=o + ^ (-5^ r=o r)+... + ^-(r'F(r)) „ = 0. (3.7) 



The mean dissipation rate of the potential enstrophy is the leading order term of the 
time derivative: 

|(a(0)) = |^)=2|(Q} = -2, Q ^ ^ ^ (3.8) 

where it is assumed that the higher-order contributions to the time derivative in Eq. I|3.6|l 
vanish as r — > in the statistically steady state. This is an alogous t o the assumption 
made to derive the Karman-Howarth equation for helicity llKurienl i20()$\ ). Next use 
l|3.7|l - l|3.8(l in l|3.6|l . multiply by r 2 throughout, and integrate with respect to r to find 

F{r)~-~e Q r, r -» 0. (3.9) 

The constant of integration is zero assuming that the scalar function F(r) remains regular 
as r — > 0. Alternatively, we can write the '2/3-law' for the third-order correlation of q 
and velocity as: 

2 

(q(x)q(x + r)(u(x) - u(x + r)) • r)) = -^ e Q r (3.10) 

where eq is obtained from Eq. (|2.3|l in the large Ro, large Fr limit, 

EQ = ~{v{q V 2 w • Vp) + K(g V 2 Vp • (3.11) 

The equation (|3.10() is the potential enstrophy counterpart to the scaling laws presented 
in Eq. 1)1. l|l . As for the helicity case, cq is not sign definite. As we have discussed above, 
this fact in itself does not exclude the possibility of an inertial range. 

3.2. The quasi- geostrophic limit 

The quasi-geostrophic (QG) limit is attained in the high-rotation, large-stratification 
(low Ro, low Fr) limit corresponding to 

N 

Ro = e—, Fr = e, Re > 0(1), Pr fixed as e -> 0. (3.12) 
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Case Ro Fr q viscous-diffusion terms 



IN 1 



w ■ Vp fie" 1 ^ {{qp'&M - q'pd]uJi)) 



c 



+(RePr)~ 1 dr l {(qu'ifyp - ?Wf p)) 
6 O(l) |^ 2( J R e Pr)- 1 9 r 2 j {qg') 



iv 0(1) e ^ 3 2Re- 1 d? j (qq'} 

v i 0(1) Fr-^s+Widip Re-^dqp'dfa'i-q'pdjUi)) 

+(RePr)- 1 d ri ((q< J j' t d],p' - gW|p)} 

vi 0(1) i ftj-'l+Uiftp Re- 1 dr i {(qp'd?,w' i -q'pd]u H )) 

+ {RePr)- 1 d r% {{quj' az d 2 jl p' - q'uja^p)) 
Table 1. The form of q and the viscous-diffusion terms of 12.91 in various cases. 



After change of variables, ensemble averaging and using homogeneity, the leading order 
0(l/e 2 ) terms from J22J are 

9 J±fL - d n (q q '{u % - «{)> = jfe-^a*^) - <4> + 

J* ^_ flr<< __>__ { ^ (wjl _> + ^(^_>}j, «~«»--^. (3.13) 

In this limit, the viscous-diffusion terms on the right-hand side of (|3.13() are all second- 
order derivatives with respect to r. Therefore they can be considered localized at small 
scales, allowing an inertial range in which viscous-diffusion contributions are separated in 
scale from transfer contributions. To focus on QG scales rather than small scales , we scale 
out the dependence on / and N by using a stretched ^-coordinate, z* = (N / f)z llCharnevI 
(119711) ). Then the scaled q is given by q = w 3 — dp/dz*, and all z-derivatives (f/N)d/d z 
in i|3.13|l become d/dz*. Such vertical stretching removes anisotropy in the QG scales; 
one may imagine pancake-shaped eddies, stretched to become spherical eddies. Following 
sections 13.1.11 and 13.1.21 one may next derive the Karman-Howarth equation analogous 
to Eq. (|3.6|l . and then a scaling law for the third-order mixed correlation function in an 
isotropic range of scales, identical to Eq. (|3.10() . However, the separation distance r is 
measured in the stretched coordinates, and the mean dissipation rate for QG is given by 



dx 2 dy 2 f 2 



4. Summary and Discussion 

Our results for various limits of Ro and Fr are summarized in Tabled In cases (ii-iv), 
an inertial range of scales dominated by flux of potential enstrophy may be realized at 
finite Re because the viscous-diffusion terms are confined to small scales. For QG flow 
with Ro,Fr — ► (case ii) and Re —> oo, the 2/3-law describes locally isotropic flow 
in a coordinate system with vertical stretching. In the remaining three cases (i, v and 
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vi), the possibility of a potcntial-cnstrophy inertial range at finite Re is not as clear 
because the flux and viscous-diffusion contributions can in principle intermingle at all 
scales. Nevertheless, the 2/3-law Eq. (|3.1U|> is obtained for the Ro, Fr — > oo case (i), 
again for locally isotropic flow in the Re — > oo limit. In the derivation of the 2/3-law 
for both cases (i) and (ii), there is a prescribed sequence of limits: Re — > oo with Pr 
fixed, followed by the separation r — > 0. In case (i), the 2/3-law is additional statistical 
information describing 3D turbulent flow with a passive scalar p. A point of difference 
from Eq. (1.1b) is that the 2/3-law for q contains information about the geometry and 
structures of the flow; recall that u> evolves as a line el ement and V p evolves as a surface 
element, and thus that q evolves as a volume element l|Ertell l|l942() '). 

The 2/3-laws are derived assuming local isotropy for some range of scales: "sufficiently 
small" scales for nearly isotropic flow with Ro, Fr — > oo; larger QG scales for which 
vertical stretching removes anisotropy in the case of Ro, Fr — > 0. It is important to note 
that even if local isotropy is not strictly reali zed, there is a l ways a n isotropic component 
of the flow that could obey the 2/3-law (see iTavlor et all l|2003j) '). We expect the 2/3- 
laws to become part of the theoretical foundation of rotating and stratified flows, and 
a benchmark for high-Reynolds number physical and numerical experiments, just as 
Eqs. I|l.l|) provide a foundation for isotropic turbulence research. 

An eventual goal is to elucidate the connectio n between energy and potential enstrophy 
in different parameter limits, as pioneered by ICharnevI l|1971|) for QG dynamics. It is 
instructive to write the conservation laws for energy and potential enstrophy: 

|<|v| 2 + p 2 > = i?e- 1 (V 2 |v| 2 ) + (fePr)- 1 (V 2 p 2 ), 
= ife-^f r-^gVVj) - (iVp ■ V 2 u;» - Pr' 1 {Ro~ l (qV 2 ^) + (q u, ■ V 2 Vp))), 

where q = oj ■ V/5 + Ro~ 1 dp/dz — Fr~ 1 LU3. The energy conservation law does not depend 
on either Ro or Fr, while the conservation law for potential enstrophy does (see also 
Table Thus potential enstrophy may constrain energy transfer among different scales 
in cases other than the QG limit, an exciting possibility for future research. 
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